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SUMMARY 

The transfer of heat from a metallic ribbon to water, under tran- 
sient conditions, was studied for initial water temperatures from 76° to 
203° F and average heat- generation rates per square foot of surface area 

of 3, 10, and 20x10^ Btu per hour per square foot. The tests were con- 
ducted at atmospheric pressure under pool boiling conditions. The power 
surge duration was held constant at 30 milliseconds. 

An analysis assuming time periods of such short duration that free- 
convection effects are negligible was successful in predicting ribbon 
temperatures as a function of time and heat- generation rate in the non- 
boiling region. A simplified analysis assuming the ribbon to be com- 
pletely insulated was found to be applicable during time periods in the 
nonboiling region of short duration. In addition, the analysis was 
found to be applicable for the entire surge duration under conditions of 
heat-generation rates well above those required for film boiling. 

Transient critical heat flux values up to 3XL0 6 Btu per hour per 
square foot were obtained for water at atmospheric pressure. The ob- 
tained data indicated an increase in critical heat flux with increase in 
the degree of subcooling for a constant heat- generation rate, and an 
increase in critical heat flux with increase in the rate of heat genera- 
tion for a constant degree of subcooling. 


INTRODUCTION 

Of current interest in the field of nuclear reactor design is the 
problem of fuel element burnout due to a sudden large power surge . In 
the advent of an exponential power surge resulting from a step change 
in reactivity, where a subcooled liquid is used as a coolant, the wall 
temperature rises rapidly until it has exceeded the saturation tempera- 
ture of the fluid. As a consequence of the temperature excess, surface 


2 


NACA TN 3885 


(nucleate) boiling occurs resulting in an increase in the rate of heat 
transfer. Provided that the heat- generation rate is still greater than 
the rate at which heat is being transferred, vapor formation increases 
until a point is reached at which the bubbles coalesce (film boiling) . 

The increased thermal resistance, due to the film of vapor, causes a 
rapid decrease in the heat- transfer rate and a corresponding rise in 
wall temperature, resulting in the eventual melting of the fuel element. 

To determine the minimum safe reactor period (the minimum time taken 
for the power to increase by a factor of e without endangering the 
reactor) which is allowable with a given control system design, a knowl- 
edge of the heat-transfer characteristics during the surge duration is 
required (ref. l) . The purpose of this paper is to give an insight into 
the transient heat- transfer characteristics of a liquid- cooled system 
during a sudden power surge. The experimental system consisted of an 
electrically heated nickel ribbon immersed horizontally in a pool of 
water, with a sudden large power surge created by a step change in 
voltage. 


APPARATUS AND PROCEDURE 
Main Circuit 

The main circuit (fig. l) contained a nickel ribbon suspended hori- 
zontally by means of two copper rods in a pool of boiled distilled water, 
a calibrated resistance, and a bank of four 12-volt storage batteries. 

The nickel ribbon, which was 3 inches by 0.15 inch by 0.005 inch, was 
heated electrically by a current pulse originating in the storage bat- 
teries. The two copper rods acted as heat sinks and remained at the 
same temperature as the mass of water in which the ribbon was submerged. 
By means of a dual- channel oscillograph and an oscillograph-record 
camera, the voltage drop across the copper rods, which was essentially 
that across the ribbon, and the voltage drop across the shunt were meas- 
ured and recorded as a function of time. Instantaneous values of average 
temperature and heat- generation rate were obtained from these measure- 
ments since the ribbons had previously been calibrated as resistance 
thermometers. Variations in liquid subcooling and heat- generation rate 
were accomplished as follows. For a given series of runs, the voltage 
drop across the ribbon was held constant and the degree of subcooling 
varied. The voltage drop was then increased and held constant while the 
degree of subcooling was again varied. This procedure was followed until 
a voltage drop setting was reached at which the ribbon failed. For a 
given voltage drop setting, the initial heat- generation rate decreased 
by a maximum of 13 percent as the degree of subcooling decreased. Also 
for a given run, the heat-generation rate decreased by as much as 50 per- 
cent with time. In order to distinguish between the various levels cor- 
responding to a given voltage drop, an average value for each is used 
throughout the report. 
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Timing Circuit 

The timing circuit (fig. 2) regulated the following items hy means 
of a group of relays and a motor-driven hank of timing wheels: a 16- 

^•^■Tlimeter camera, a flash bulb, the driven sweep of the oscillograph, 
and the power surge duration. The camera, operating at a rate of 5000 
frames per second, recorded the boiling action on the ribbon. The indi- 
vidual frames on the film were related to the voltage traces on the 
oscillograph by means of 1000-cycle-per-second timing marks applied 
simultaneously to the z-axis of the oscillograph and an argon timing 
lamp within the camera. 


Methods of Calculation 

The basic data needed for the determination of the system heat- 
transfer characteristics are the simultaneous voltage drops across the 
ribbon and a calibrated resistance in series with it. The voltage drops 
were determined as a function of time by means of a dual- channel oscil- 
lograph and an oscillograph-record camera. Calibration was obtained on 
each print by means of a known square wave voltage for each beam, obtain- 
able from the oscillograph. 

Rate of heat generation . - The instantaneous rate of heat generation 
was determined by the product of the voltage drop across the ribbon and 
the line current which is obtainable from the voltage drop across the 
calibrated resistance 


P = 3.413 vl (l) 
Symbols are defined in appendix A. 

Average ribbon temperature . - The ribbon resistance was determined 
from voltage and current measurements by the assumption of a pure resis- 
tance circuit 


R(0 = j 

Calibration of the nickel ribbon as a resistance thermometer determines 
the instantaneous average ribbon temperature. The calibration curves of 
resistance as a function of temperature are given in figure 3. 

Rate of heat transfer . - The rate of heat generation is given by 
equation (1) . The rate of heat absorption is given by 

d0 w 

wc p at" 
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Neglecting end losses, the difference between these two quantities must 
equal the instantaneous average heat flux 


q 


1 

S 


^3.413 vl 


WC 


d e 

P dt 


w 


( 2 ) 


with d0 w /dt obtained from a plot of average ribbon temperature against 
time. The Cp term was obtained as a function of average ribbon tem- 
perature by the relations (ref. 2) 


c p — 4.26 + 0.00640 T (T = 273° to 626° K) 
c p = 5.99 + 0.000905 T (T = 626° to 1725° K) 

(cal/(deg) (mole) 

To justify the neglect of end losses in equation (2), an analytical 
expression is developed in appendix B to illustrate the type of temper- 
ature distribution obtainable along an internally heated slab of finite 
length, each end of which is in contact with a heat sink, if the mode of 
heat transfer is assumed to be solely by conduction along the length of 
the slab. The length of the slab was taken to be 3 inches, and a heat- 
generation rate comparable to those resulting experimentally was used. 
From the results shown in figure 4, the heat loss through the ribbon 
terminals was calculated to be 0.26, 0.28, and 0.3 percent at 2, 5, and 
8 milliseconds, respectively. 


Heat- transfer coefficient . - The heat-transfer coefficient, based 
on average ribbon temperature rather than surface temperature, is de- 
fined by the equation 


q = h(0 w - 9 0 ) 

The instantaneous average heat- transfer coefficient was therefore deter- 
mined by the value of q previously obtained and the corresponding 
temperature difference. 


RESULTS AND DISCUSSION 

Because of the very short time periods employed, it might be ex- 
pected that average ribbon temperatures could be predicted with a fair 
degree of accuracy in the initial moments of the period by assuming the 
ribbon to be completely insulated. Figure 5 illustrates the divergence 
of the experimental data from that assumption. The maximum deviation of 
the average ribbon temperature in the nonboiling region was 43 percent. 
In the nucleate boiling region, as the initial system temperature was 
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increased, the deviation from the assumption of complete insulation in- 
creased. In the partial and complete film boiling regions, the deviation 
decreased with increasing initial system temperature, approaching a value 
of constant magnitude. In appendix C, an analytical expression for the 
average ribbon temperature is developed assuming the ribbon to be com- 
pletely insulated. 

Because of the fair agreement found in the nonboiling region, the 
possibility arises that during the short time period used free- convection 
effects are negligible and, therefore, that heat transfer to the sur- 
rounding medium is by the process of conduction. In appendix B, it is 
shown that heat loss through the ribbon terminals was negligible. Since 
the ratio of ribbon width to ribbon thickness was 30, it may be assumed 
that the major portion of the heat flow was through the face of the rib- 
bon. Thus, the problem may be considered one-dimensional. In appendix 
D, analytical expressions are derived for the temperature distribution 
through the ribbon thickness and into an infinite medium. The resulting 
distribution, using the same heat- generation rate as in figure 4, is 
given in figure 6. Although it may be seen from figure 6 that the tem- 
perature distribution within the ribbon was quite flat, for completeness, 
an analytical expression is derived in appendix E for the average ribbon 
temperature . 

In figure 7, the average ribbon temperatures as derived in appendixes 
C and E are compared with the experimental data. The agreement with 
appendix E is excellent, bearing out the validity of the assumption. 

Such agreement also indicates that in the nonboiling region, the average 
temperature of the ribbon did not deviate greatly from its surface tem- 
perature. Thus, in the nonboiling region, heat- transfer coefficients 
would correspond to those based on surface temperature measurements. In 
the nucleate boiling region, there is some uncertainty as to the devia- 
tion of surface temperatures from average temperatures. However, since 
the average temperatures measured seem to correspond in magnitude to 
surface temperatures normally encountered in nucleate boiling, the cal- 
culated heat-transfer coefficients should be at least usable. Figures 
8 and 9 are typical records of the variation of heat flux and heat- 
transfer coefficient with time and average temperature difference, cor- 
responding to the temperature- time curves in figure 7. 


Nucleate Boiling Region 

Average ribbon temperature . - No reliable prediction of the ribbon 
temperature could be made in the nucleate boiling region. The ribbon 
temperatures ranged from 220° to 280° F depending upon the rate of heat 
generation and degree of subcooling. The average of the ribbon temper- 
atures in the nucleate boiling region was 250° F from which the deviation 
was ±16 percent. 
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The transient "behavior of the average ribbon temperature in the 
neighborhood of the inception of boiling is such that the temperature 
gradient decreased, deviating from the analysis of appendix D. Depending 
upon the rate of heat generation, degree of subcooling, surface condi- 
tions, and gas content, a peak in the temperature curve may result, such 
that a discontinuity is obtained in the temperature gradient. Figure 10 
illustrates this condition. The discontinuity has also been observed 
in an investigation similar to the present one by using an exponentially 
increasing heat- generation rate (ref. 3). 

Heat transfer . - The inception of nucleation disturbs the thermal 
boundary layer, and the resultant turbulence causes a sharp increase in 
the rate of heat transfer. Depending upon the rate of heat generation 
and the degree of subcooling, the bubbles may expand and collapse, leave 
the ribbon entirely, or coalesce to form a film of steam completely sur- 
rounding the ribbon. Under the first two conditions, the rate of heat 
transfer approaches some steady- state value. 

Excerpts from the film record of the boiling action corresponding 
to figures 7(c), 8(c), and 9(c) are shown in figure 11. Comparison of 
the variation in the heat-transfer characteristics with boiling action 
may be made, as the elapsed time between the initiation of the surge and 
the exposure of the particular frame is indicated. 

By reference to figure 8(c), it may be inferred that surface boiling 
began at approximately 1.6 milliseconds after the initiation of the power 
surge, since at this time its effects on the heat-transfer coefficient 
were noticeable. Thus, the first frame shown in figure 11 was exposed 
approximately 1 millisecond after the inception of boiling. The bubbles 
which seem to be in the process of leaving the ribbon were located on 
the walls of the pyrex beaker and should not be associated with the vapor 
formation on the ribbon. Figure 7(c) shows that in the time period from 
2 to 3 milliseconds, the temperature curve deviated from that predicted 
in appendix E. Between 2.7 and 3.1 milliseconds, the bubbles grew radi- 
ally retaining their hemispherical shape. As they grew and as new bub- 
bles formed, more surface area was covered by vapor. Thus, at some 
point the thermal resistance of the partial vapor film overcame the de- 
crease in thermal resistance resulting from the turbulence caused by the 
growth of new bubbles and expansion of old bubbles, with the result that 
at 3 milliseconds the rate of heat transfer began to decrease (fig. 8(c)). 
From figure 7(c) it may be seen that an inflection occurred at the same 
instant, and the temperature gradient began to increase. The formation 
of new bubbles and growth of older ones continued with the result that 
the bubbles partially coalesced, giving rise to what is known as the 
partial film boiling region. At 6.7 milliseconds, the ribbon seemed to 
be surrounded completely by bubbles. However, large individual bubbles 
in the outer regions were still discernible. At 8.4 milliseconds, the 
bubbles had coalesced to the extent that a film of vapor surrounded the 
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ribbon. Figures 7(c) , 8(c) , and 9(c) show that the rate of heat transfer 
leveled out and began to increase and the temperature gradient slowly 
decreased. These characteristics are associated with what is known as 
the stable film boiling region. The absence of data after a time dura- 
tion of 18 milliseconds was due to a partial obscuring of the voltage 
trace at that time plus the fact that the ribbons were not calibrated 
past 800° F. After 18 milliseconds the heat- generation rate continued 
to level out, and it may be assumed that the ribbon temperature contin- 
ued to rise. The actual length of the surge duration was 26 milliseconds. 
Referring again to figure 11, it should be noticed that there was only 
a very thin film covering the bottom of the ribbon; thus the expansion 
of the film occurred on the upper face. At 13.0 milliseconds, there was 
a noticeable disturbance of the film close to the upper face of the rib- 
bon. This disturbance grew until at 25 milliseconds, the entire film 
seemed to be tearing away. An interesting feature of the series of 
frames from 25 to 35.4 milliseconds was the seemingly thin membranelike 
sheet of film which connected the leaving vapor with the thin film re- 
maining on the ribbon. At 45.3 milliseconds, the vapor bubble had com- 
pletely dissociated itself, leaving thin films of seemingly equal thick- 
ness on the top and bottom of the ribbon. 

In general, the first observable density disturbance occurred at the 
moment at which the ribbon temperature began to deviate from its pre- 
dicted value. The disturbances (vapor bubbles), when first observed, 
were hemispherical in shape and had a height of the order of 0.005 inch. 
The disturbances expanded rapidly and broke up the thermal boundary 
layer, the resulting increase in heat transfer causing the slope of the 
temperature- time curve to approach zero. 

The subsequent expansion of the bubbles perpendicular to the ribbon 
was dependent upon the degree of subcooling of the surrounding medium. 

For low degrees of subcooling, the bubbles retained their hemispherical 
shape while expanding. As the degree of subcooling increased, the rate 
of expansion perpendicular to the ribbon decreased and the bubbles tended 
to grow horizontally along the length of the ribbon. In both instances, 
many of the individuals combined to form larger bubbles which instead of 
collapsing completely, as in the case of single disturbances, reacted to 
present a pulsating type of motion. Although these bubbles did coalesce 
to some extent, the positions of the individuals were still discernible. 

It might be remarked that the growth of many of the disturbances 
perpendicular to the ribbon was essentially complete at the time they 
were first observed. Since the film rate was 5000 frames per second, 
the bubbles had grown to vertical maturity in less than l/5 of a 
millisecond. 
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No correlations were obtained in the nucleate boiling region for 
either heat- transfer coefficients or rates of heat transfer. In si l 
cases, the inception of boiling was marked by a sharp increase in the 
heat- transfer coefficient. The behavior of the coefficient under tran- 
sient conditions is interesting. Although considerable error could have 
been introduced by the necessity of having to measure temperature gra- 
dients, the results in general appear to be quite consistent. As noted 
previously, h should decrease as the thermal boundary layer builds up. 
Generally, such was the case. Under these conditions, the coefficient 
decreased until the boundary layer was disturbed by the formation of 
small vapor bubbles. The resultant sharp increase in h, caused by the 
turbulence, gave rise to a minimum condition. In most cases this minimum 
was so well defined, especially for a plot of h against G v - Q q , that 
the corresponding average ribbon temperature and time at which nucleation 
occurred may be obtained directly, without photographic equipment. 

Values of nucleate boiling heat- transfer characteristics obtained 
in this investigation are tabulated in table I. Only those values which 
occur prior to the critical heat flux are listed for cases in which film 
boiling dominates. For cases in which a negative temperature gradient 
resulted after the inception of boiling, the values are listed only for 
positive or zero gradients. 


Film Boiling Region 

Average ribbon temperature . - No attempt was made to predict the 
ribbon temperatures in the film boiling region, although it appears pos- 
sible that an analysis of the same type as that employed in the nonboil- 
ing region is possible. The problem would be that of heat transfer by 
conduction, from the ribbon through a film of steam and into an infinite 
medium. 

Upon transition from the nucleate boiling region to the film boiling 
region, the ribbon temperature rose rapidly and if heating were continued 
the ribbon would have been destroyed. As shown in figure 5(c), as the 
degree of subcooling was decreased, the rate of heat absorption approached 
the rate of heat generation. Thus, at high heat-generation rates and 
fluid temperatures near saturation, the ribbon temperature corresponded 
closely to that which would have occurred had the ribbon been completely 
insulated. 

Heat transfer . - In the transition from nucleate to film boiling, 
the vapor bubbles expanded as previously described until they coalesced. 
The ribbon was completely surrounded by a number of vapor bubbles which 
had merged together, although the large individual bubbles were still 
discernable. Because of the relatively thick layer of bubbles, the rate 
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of" heat transfer was decreased and the ribbon temperature consequently 
rose. As the temperature continued to rise, the film increased in thick- 
ness until a condition of such instability was reached that the film 
exploded, shooting off bubbles of enormous size relative to the ribbon. 

A very thin pulsating film in which individual bubbles were no longer 
discernible immediately formed and the transition to stable film boiling 
was complete . 


Critical Heat Flux 

In the region of transition from nucleate to film boiling, there 
is a maximum rate of heat transfer. This value is usually termed the 
critical or maximum heat flux. 

Experimentally obtained transient critical heat flux values are 
tabulated in table II. The third to ninth values indicate that for con- 
stant heat- generation rate the critical heat flux decreased as the degree 
of subcooling decreased. The second and eighth or first and ninth 
values indicate that, for a given degree of subcooling, the critical heat 
flux increased as the heat- generation rate increased. 

The bubble formation in the neighborhood of the critical heat flux 
corresponding to figures 7(c), 8(c), and 9(c) is illustrated in figure 
11. From figure 8(c) it may be seen that q c occurred at 3 milli- 
seconds. Between 2.9 and 3.1 milliseconds the bubble formation expanded 
and covered a larger area of the ribbon. The increased coverage of rib- 
bon area by bubble formation offset the turbulence created by the expan- 
sion of old bubbles and the formation of new bubbles with the result 
that the thermal resistance increased and consequently the rate of heat 
transfer decreased. 


SUMMARY OF RESULTS 

The following results were obtained from the experimental and 
analytical investigation of transient water pool boiling resulting from 
a sudden large power surge: 

1. For time periods in the nonboiling region up to 30 milliseconds 
(the maximum time duration investigated), free-convection effects were 
negligible and heat transfer therefore occurred by the process of con- 
duction. Thus, ribbon temperatures may be predicted as a function of 
time and rate of heat generation up to the inception of boiling. 

2. For time periods in the nonboiling region of less than 14 milli- 
seconds, the maximum deviation of the average ribbon temperature from 
that obtained by assuming the ribbon to be completely insulated was 43 
percent. 
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3. In the nucleate boiling region, the average ribbon temperature 

was 250° F. The deviation from this average was ±16 percent. 

4. For heat- generation rates much greater than those required for 
film boiling, the ribbon temperature was predicted over the entire range 
(nonboiling and nucleate and film boiling) by the assumption that the 
ribbon is completely insulated. 

5. Transient critical heat flux values up to 3x10^ Btu per hour per 
square foot were obtained for water at atmospheric pressure. 

a. For a constant heat- generation rate, the critical heat flux 
increased as the degree of subcooling increased. 

b. For a constant degree of subcooling, the critical heat flux 
increased as the rate of heat generation was increased. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, Sept. 17, 1956 
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APPENDIX A 
NOMENCLATURE 

constants 

constants 

value of ordinate on plot of heat -generation rate per unit 
volume- time, Btu/(hr)(cu ft) 

constants 

specific heat, Btu/(lb)(°F) 
heat source, Btu/(hr)(cu ft) 
error function 

complementary error function, 1 - erf 

heat- transfer coefficient, Btu/(hr)(sq ft) (°F) 

line current, amp 

thermal conductivity, Btu/((hr)(sq ft)(°F)/ft)) 
length of ribbon, ft 

slope on heat generation rate per unit volume-time plot, 
Btu/(hr B )(cu ft) 

integer, 1, 2, 3, . . . 

power or rate of heat generation, Btu/hr 

thermal energy generated in a system, Btu 

thermal energy transferred to a system per unit time per 
square foot of surface area, Btu/(hr)(sq ft) 

critical heat flux, Btu/(hr)(sq ft) 

ribbon resistance, ohms 
m + sb 
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surface area of ribbon, sq ft 

variable introduced by application of Laplace transformation. 
If f(t) is a known function of t for values of t>0, 
its Laplace transform f(s) is defined by the equation 

f(s) = r e" st f(t) dt (ref. 5) 

Jo 


absolute temperature, °K 
time, hr or millisec as specified 
volume of ribbon, cu ft 
voltage drop, volts 
mass of ribbon, lb 
coordinate system, ft 
thermal diffusivity 

k/aVz 

h/(f>z + l>l) 

(32 " 3i)/( 32 + 3i) 

temperature, °F 
initial system temperature, °F 

coolant saturation temperature minus average coolant temper- 
ature, °F 

average ribbon temperature, °F 
Laplace transform of 0 



density, lb/cu ft 

l/2 thickness of ribbon, ft 


— , sq ft/hr 
P 
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Subscripts : 

1 region of ribbon 

2 region of medium surrounding ribbon 


* 
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APPENDIX B 


TEMPERATURE DISTRIBUTION ALONG LENGTH OF RIBBON 


This analysis is intended to illustrate the type of temperature 
distribution (fig. 4) ohtainahle along an internally heated finite length 
of ribbon suspended horizontally by means of two copper rods, if the 
mode of heat transfer is assumed to be solely conduction along its length. 
The model may be described as follows: 

(1) It is an infinite slab of finite length 2, and each end is in 
contact with a heat sink. 

(2) The slab has conductivity k, density p, specific heat Cp, 
and a heat source d. 

(3) The temperature distribution satisfies the following equations: 


The experimental data indicate that d may be expressed as a linear 
function of time. Equation (Bl) thus becomes 


d 2 fl(x,t) _ P c p d9(x,t) _ _ d 


dx 2 k ~ ^ k 


(Bl) 




d 2 6>(x,t) P c p d6 (x, t) _ (mt + b) 
^ x 2 k at ~ k 


(B3) 


Applying the Laplace transformation to equations (B2) and (B3) 


d2 S( x ' s ) - n^xjs) = -r 


dx' 


(B4) 


where 




and r = o 

ks 2 



(B5) 
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Making use of equations (B2), equation (B5) becomes 


0(x,s) = — £ 


4 


-[IX 




1 + e“^ Z 1 + e^ l j 
Equation (B6) may be written as (ref. 4) 

0(x,s) = 


(B6) 


,jl - |V“ X + e-^< l - x >][l + 

Expanding the quantity ^1 + e yields 

(l + e-^y 1 = 1 - e-^ + e" 2|il - ... + (-l) n e^ nZ = (-l) n 

rt=0 


where 0,1,2,... and e < 1 

Substituting equation (B8) in equation (B7) results in 


(B7) 


-iml 


(B8) 


0(x,s) « 


- [■ 


•^ x + e "H(7- x ) 


>]£ M> n *'4 

J n=0 J 


(B9) 


Expanding equation (B9) and substituting for r 

°° 00 

0(x,s) - + b)< 1 - (-l) n e~^ nZ+x ) - yy (-l) n e ”^C Z ( n+1 ) -x J j 


n=0 


n==0 


(BIO) 


Taking the inverse transform and simplifying yields 


0(x,t) = ~ (mt + 2b) - 

16omt^ i 4 erfc x 


x +) (-l) n+1 i 4 erfc 21 ~ x - i 4 erfc + x l 

2-y/at Z ^ L 2-y/ at 2-y/ at J J 


iTOi 2 erfc — ^ 


2-v/at £ — ' 

v n=l u 


i 2 erfc ~ — -i 2 erfc nl 


2-y/at 


nl + x I 

2 A /at J 


(Bll) 
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APPENDIX C 


SIMPLIFIED DERIVATION OF RIBBON TEMPERATURE 


It is assumed that during short time periods, heat transfer from 
the ribbon can be neglected in the nonboiling region and, therefore, that 
the rate of heat absorption equals the rate of heat generation; the rib- 
bon is completely insulated. 


a(Q/v) „ 

dt 


d G 


pc 


w 


P dt 


(Cl) 


Equation (Cl) may be written in integral form as 



(C2) 


The ribbon temperature may be evaluated from equation (C2) by determining 
the area under the curve resulting from a plot of d(Q/V)/dt against t. 

The heat- generation rates obtained experimentally from a step change 
in voltage were found to be representable as a linear function of time, 
that is, 


a (aA) . mt + b 

dt 


(C3) 


Substituting equation (C3) into equation (C2) gives 

r e v ! r b 

d0 w “ 77" / ( mt + b ) dt 

Je 0 pc pJo 


(C4) 


Assuming p and Cp constant, equation (04) becomes upon integrating 

t 


= 0 o + 2pc t 


(mt + 2b) 


(C5) 
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APPENDIX D 


TEMPERATURE DISTRIBUTION THROUGH RIBBON INTO AN INFINITE MEDIUM 


This analysis is intended to illustrate the type of temperature dis 
tributions (figs. 5(a) and (b) ) obtainable along the finite thickness 
of an internally heated ribbon in contact with an infinite medium, if 
the mode of heat transfer is assumed to be conduction through the ribbon 
and into the surrounding medium. The ribbon may be described as follows 

1. It is an infinite slab of finite thickness in contact with an 
infinite medium. 


2. The slab (region l) has conductivity kp, density pp, specific 
heat Cp and a heat source d. 


3. The infinite medium (region 2) has corresponding constants kg, 

Pg, and c p,2 no heat source. 

4. The temperature distributions, 0p(y,t) and 02(y>t), satisfy the 
following equations: 


S 2 0p(y,t) 

dy 2 


Pl c p,l d0i(y,t) 
kp St 


_d_ 

k l 


0 < y < t 


d 2 0 2 (y;t) P2 c p,2 

Sy 2 ~*2 5t 


0p(T,t) 

= 0 


r 

\ ( 

ae 2 \ > 

V 1 ^y , 


dy ) 


'y** J 



(Dl) 


(D2) 


(D3) 


(D4) 


y 
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The experimental data indicate that d may he expressed as a 
linear function of time. Equation (Dl) thus becomes 

d 2 0, (y,t) p, c . de (y,t) , # , 

- B?— ^ = - - — (mt + h) 0 < y < f (D5) 

a y 2 k i s k i 

Applying the Laplace transformation to equations (D2), (D3)> (IH), 
and (D5) results in 


ci 2 0 n (y,s) 


2 - 4 9 l (y ' S) " _r 


dy 

d*e 2 (y,s) 

O ^2 9 2^ y ’ S ^ “ 0 




(D6) 


(D7) 


where 


1.2 a m + sb 

, and r = — 5- 

k^s^ 


- -S— . an o = ^ 

h,2 <^,2 1,2 P 1j2 c p ,1,2 

The solutions of equations (D6) and (D7) may be written 

5(y,s) - A 1 (s)e"^ y + B 1 (s)e' 1 l y + ~ 

0 2 (y,s) - A 2 (s)e" M ' 2y + B 2 (s)e^ 2y 
By use of equations (D4), equations (D8) and (D9) become 


9l(y^s) = ApCs)^! 7 + e ^ iy ) + 


r 

3 


9 2 (y,s) = ^(sje 


-P 2 y 


(D8) 

(D9) 

(DIO) 

(Dll) 


Making use of the boundary conditions in equations (D3), the constants 
A p and Ag in equations (DIO) and (Dll) are found to be 


Ai = 


rk 2 n 2 


|^p. 2 k 2 (e^ lT + e ^ lT ) + n 1 k 1 (e 111 - e )] (D12) 


rk 

A 2 ** ^ 


( e^ - 




1 V e- - e 1 ') j-^ k2Ce Ult + + . e -^)] 


-1 


-P2 t 


(D13) 
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Substituting equations (D12) and (D13) into equations (DIO) and (Dll) 
results in 


0i(y,s) ^ ~Ji . k 2 u 2 


^i y + e -^i y 


Ml | [^(e^ + e"Ml T ) + mkiCe^ - e“Ml^) 


rk-, / -M 2 y\ 

e 2 (y ’ s) - crlz^J 


(D14) 


e Ml*T _ e -HlT 


+ e Ml T ) + f^k]_ (e^ lT - e”Ml T ) 

Simplifying equations (D14) and (D15) gives the equations 

Sl(y,s) - ^|l - 3^ [e-"l (T - y ) + e-“l(-' +y )](l ♦ {e'H 


(D15) 


V 

*»<*>•> - - *' 2tIlT )(^ * c .** 1 ’)' 1 


(D16) 

(D17) 


where 


k l,2 P 2 ' P 1 

h,2 « 375- 5 * P 2 + Si 

Expanding the quantity (l + £e -2 ^-*- T ) yields 

(l + te- 2 ^)' 1 . 1 - ? e - 2 W T + t V 4 M- T - ... + (.1)“ ^-Z^n 

m 

Y, (-l) n £ n e" 2 W- Tn 


n^O 


(D18) 


for n«0,l>2,... and £ 2 e ' 4 Ml T < 1 
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Equations (D16) and (D17) then become 

5 l(y , B ) . . r [e-^-y) + Y, (-1)“ 5 n e- 2 “l TO j 




(D19) 


Q?(y>s) - 


rk]_a : 


1/2 


bn E 


s(P 2 + Pi) 


e -|i 2 (y-T) _ e -|i 2 (y-T)-2mT 




(-i) n c 


n /'n e _ 2M-i T n\ 


r^O 


(D20) 


where r = p 2 + Pi 

Expanding, taking the inverse transform, and simplifying equations (D19) 
and (D20) result in 
a^t 

®l(y>t) = 2 v (mt + 2b) - 


V 1 
16a- 


4a^bt 


1/2 


2^Llly ( .i)n c n(i4 erfc + l4erf<; h(2yl),y^ 

kl \X l 2 -Va^t jj 

/X 


(-l) n ^i 2 erfc 


[r(2n+l)-yj 2 [x(2n+l)+y] | 


+ i erfc 


e 2 (y-t) - p/m h E ( ‘ i)n?n( 4 t) i erfo 

V n=0 v. 

[_a}/ 2 (y - t) + 2a 2 / 2 x(n + l)j| 


2- v / a i't ' " """" 2-y/a^t 

r 1 / 2 , . i/s 1 

L a l (y - *0 + 2a 2 TnJ 


(D21) 


2 - \Z a 2 a l’ t 


i^erfc 


2-y/ a 2 a it 


+ 


b (-l) n ?; n (4t)^ i 2 erfc 
n=0 


jjx^ 2 (y - t) + 2a^^ 2 -rnj 


i^erfc 


2-y/aga^t 

[a^ 2 (y - t) + 2a 2 ^ 2 T(n + l)J 
2-y/a^a^t 


(D22) 
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APPENDIX E 


AVERAGE RIBBON TEMPERATURE 

The equations developed in appendix D for the temperature distri- 
butions may be integrated as follows to yield an equation representing 
the average ribbon temperature: 




0l(y>s) dy 


(El) 


Substituting equation (D19) into equation (El) gives 

e v = It F% - r O l(T ' y) + e- ,1 i (T+y) ] V (- 1 ) V a- 2 “i' tn 

Un feo 


Integrating and simplifying equation (E2) yield 


dy 


ra i i r 

“ sx V " M-1 


(-l) n ^ n ^e" 2MirrT - e' 2|i l T ( n+1 )) > 

fco J 


(E2) 


(E3) 


Substituting for r, talcing the inverse transform, and simplifying 
result in _ 

y 3/2 

2nr 


a^t 


0 w “ 2k7 ( mt + 2b ) 


*“1 


V 1 * ' k-,T 


m 


^ (-l) n C n (4t) 5 / 2 i 5 erfc 
n=0 


2-y/a^t 


(-l) n 5 n («) 5 /2 l 5 erfc 

OO 

Y (-l) n ^ n (4t) 3 / 2 i 3 erfc — . 
L — ' 2-y/a^t 


_n=0 


OO 

E 

n=0 


(.l)" 5 ”(4t) 3 / 2 Perfc 

2-y/a^t 


(E4) 
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TABLE I. - NUCLEATE BOILING HEAT-TRANSFER COEFFICIENTS 


Initial 
temperature , 

op 

Average 
ribbon 
temperature 
minus initial 
system 
temperature, 

o p 

Heat-transfer 

coefficient, 

h, 

Btu 

Heat flux, 

q> 

Btu 

TErTCsq m 

Temperature 

gradient, 

d0 w /dt, 

OF/hr 

Heat- 

generation 
rate . 
d(Q/sJ 

(hr J (sq ^tJ(^F) 

dt • 
Btu 

(hr ) (sq ft) 

1 

46 

90 

2170 

1.95xl0 5 

0. 756 xlO 7 

2.92x10 s 



92 

2220 

2.04 

.684 

| 



94 

2270 

2.13 

.612 




96 

2320 

2.23 

.540 

1 



98 

2580 

2.53 

.270 

t 

166 

63 

3100 

1.95 

.810 

2.985 



66 

3590 

2.37 

.468 

2.967 


t 

70 

3600 

2.63 

.234 

2.930 

186 

51 

1530 

.78 

1.925 

3.26 



58 

3670 

2.13 

.810 

3.17 



65 

4010 

2.61 

.396 

3.12 


T 

69 

4400 

3.04 

.000 

3.04 

203 

37 

3570 

1.32 

1.640 

3.43 



43 

5720 

2.46 

.684 

3.34 


1 

47 

6620 

3.11 

.126 

3.27 


76 

149 

3950 

5.89 

3.24 

10.03 



156 

4630 

7.23 

2.16 

10.00 


1 

163 

5200 

8.48 

1.12 

9.93 

95 

129 

3880 

5.00 

4.50 

10.75 



139 

4650 

6.46 

3.24 

10.62 



144 

5320 

7.66 

2.05 

10.30 



151 

5660 

8.55 

1.30 

10.22 



153 

5960 

9.11 

.09 

10.27 

105 

117 

4000 

4.68 

4.29 

10.14 



129 

5120 

6.60 

2.63 

9.98 

i 


134 

6170 

8.28 

1.30 

9.95 

125 

97 

4210 

4.09 

5.11 

10.57 



107 

4900 

5.24 

4.10 

10.46 

\ 


120 

4980 

5.98 

3.34 

10.27 

1? 


86 

4090 

3.52 

5.61 

10.67 



98 

5800 

5.68 

3.74 

10.49 

1 

r 

106 

6900 

7.32 

2.41 

10.42 

148 

68 

3930 

2.67 

6.225 

10.57 



81 

4950 

4.01 

4.930 

10.32 



97 

6240 

6.24 

3.025 

10.14 


f 

101 

9900 

10.00 

.000 

10.00 

185 

43 

2840 

1.22 

7.31 

10.57 



65 

8720 

5.67 

7.42 

10.10 

1 


68 

11420 

7.77 

1.80 

10.00 


n 

166 

9030 

15.0 

11.10 

29.3 



184 

14120 

26.0 

2.36 

29.1 


f 

187 

15130 

28.3 

.40 

28.8 

95 

160 

6620 

10.6 

14.70 

29.6 



194 

12920 

25.1 

3.35 

29.5 



196 

14500 

28.4 

.43 

29.0 

113 

115 

11120 

12.80 

12.6 

28.9 



152 

11500 

17.47 

8.5 

28.6 

i 

r 

169 

16670 

28.20 

.0 

28.2 

131 

119 

13720 

16.32 

8.14 

26.85 

♦ 

133 

19700 

26.20 

.27 

26.55 

149 

92 

17800 

16.35 

8.64 

27.45 



112 

21000 

23.47 

2.83 

27.15 

T 


116 

21600 

25.10 

1.57 

27.15 

167 

108 

6940 

7.50 

12.88 

24.4 



138 

9520 

13.12 

8.10 

23.9 

t 


153 

10150 

15.52 

5.66 

23.1 


5 

68 

10120 

6.88 

14.05 

25.1 

f 


100 

16400 

16.40 

6.03 

24.3 

20 

3 

47 

10390 

4.88 

15.05 

24.35 



80 

11100 

8.88 

11.11 

23.50 

1 


111 

13900 

15.43 

5.52 

22.50 
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TABLE II. - CRITICAL HEAT FLUX 


Critical 
heat flux, 

Q.C > 

Btu 

Average ribbon 
temperature 
minus initial 
system 
temperature, 

“ Q o > 

°F 

Coolant 
saturation 
temperature 
minus average 
coolant 
temperature, 

^sub* 

°F 

Heat- 

generation 

rate, 

d(Q/S) 

(hr) (sq ft) 

dt * 
Btu 

(hr) (sq ft) 

3.11 X 10 5 

47 

9 

3.27 X 10 5 

7.77 

68 

27 

10.00 

28.30 

187 

135 

28.80 

28.20 

169 

99 

28.20 

26.20 

133 

81 

26.20 

25.10 

116 

63 

21 .20 

21.00 

165 

45 

23.0 

16.40 

100 

27 

24.3 

15.43 

111 

9 

22.8 
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Camer a -mo tor 
relay 


To camera 
motor 

n 

130 Volts 



24 Volts 
d.c . 



Drive -mot or 
relay 

Push button- 


To drive 
motor 



110 Volts 24 
a.c. Volts 


d.c . 



Figure 2. 


Schematic diagram 


of timing whei 



Sequence of events 

(1) Push botton in drive-motor relay circuit is 
depressed, and 16 volts are applied to drive 
motor which drives timing wheels. 

(2) Drive-motor relay closes, and push button is 
released. 

(3) Camera-motor relay closes, and 130 volts are 
applied to camera. 

(4) Flash-bulb circuit closes, and flash bulb is 
set off. 

(5) Timing-light relay closes. Oscillograph 
trace is initiated by external synchronous 
circuit, and camera timing lamp pulses. 

(6) Surge -duration relay closes, and current is 
allowed to flow in main circuit (fig. l) . 


and circuits. 
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smqo ‘souB^sfsaj uoqqfg 


Figure 3. - Resistance-temperature calibration curves for two nickel ribbons. 








Temperature 



Figure 4. - Temperature gradient along x-axis for heat -generation rate 
of (-SxiO^Jt + 8*10^ Btu per hour per cubic foot. (See appendix B.) 
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(a) Average heat -generation rate, 3xl0 5 Btu per hour per square foot. 


Figure 5. - Deviation of experimental data from that obtainable for complete insulation. 







millisec 
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12 16 
Time, millisec 


(b) Average heat -generation rate, 10x10° Btu per hour per square foot. 


Figure 5. - Continued. Deviation of experimental data from that obtainable for complete 
insulation. 
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LO 

to 




IK 



12 16 
Time, millisec 


(c) Average heat -generation rate, 20x10^ Btu per hour per square foot. 


Figure 5. - Concluded. Deviation of experimental data from that obtainable for complete 
insulation. 
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Heat-generation rate, Average ribbon temperature 

Btu/(hr)(sq ft'ixlO”^ 


MCA TN 3885 


<?• 




Time, millisec 

(a) Water temperature, 166° F. 

Figure 7. - Temperature - heat -generation-rate - time record for nickel ribbon in 
subcooled water. 



Average ribbon temperature, 




320 


280 


240 


200 


160 



CM 


120 


12 16 
Time, mi Hi sec 

(b) Water temperature, 95° F. 


Figure 7. - Continued. Temperature - heat -generation-rate - time record for nickel 
ribbon in subcooled water. 
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LO 

to 





(c) Water temperature, 203° F. 

Figure 7. - Concluded. Temperature - heat -generation-rate - 
time record for nickel ribbon in subcooled water. 
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Time, millisec 


(a) Water temperature, 166° F. 
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OJ 

cn 


Figure 8 


Heat flux and coefficient of heat transfer against time 


Heat flux, q, Btu/(hr)(sq ft) 


Heat-transfer coefficient, h, Btu/(hr)(sq ft)(°F) 
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in 



12 l 

Time, millisec 


(b) Water temperature, 95° F. 

Figure 8. - Continued. Heat flux and coefficient of heat transfer against time. 
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12 16 
Time, millisec 


(c) Water temperature, 203° F. 

Figure 8. - Concluded. Heat flux and coefficient of heat transfer against time. 
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Heat flux, q, Btu/(hr)(sq ft) 


Heat-transfer coefficient, h, Btu/(hr)(sq ft)(°F) 
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Figure 9. - Heat flux and coefficient of heat transfer against average temperature difference. 
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Figure 9. - Continued. Heat flux and coefficient of heat transfer against average temperature 
difference. 
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Figure 9. - Concluded. Heat flux and coefficient of heat transfer against average temperature 
difference. 
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Figure 10. - Temperature - time - heat-generation-rate - time record for 
nickel ribbon in subcooled water. Initial water temperature, 149 F. 
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3.9 


4.1 


4.3 



Figure 11. - Transient bubble formation resulting from sudden 
large power surge . Water temperature , 203° F . 
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32.5 



35.4 


45.3 
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C-42439 


Figure 11. - Concluded. Transient bubble formation resulting from 
sudden large power surge . Water temperature , 203° F . 
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